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ABSTRACT. 


The surface topography of a forest canopy is complex. Like the Earth’s surface, it has a 


morphology that consists of hills and valleys but punctuated by gaps. The study of the surface of the forest 
canopy has been limited, largely by access. Recent advances in remote sensing (i.e., scanning laser altim- 
etry) are beginning to provide broader views of this key interface between the atmosphere and the terrestrial 
biosphere. Our analysis of canopy topography across a 1 km? area of old-growth tropical rain forest in La 
Selva, Costa Rica—derived with laser altimetry—details patterns based on measures of depth to canopy 
surface from a given elevation and height to canopy surface above the ground. Spatial autocorrelation 
patterns of canopy height, disregarding influences of ground topography, were isotropic and significantly 
positive at scales < 50 m. The fractal dimension of canopy heights across the landscape was 1.96, indicative 
of a nearly random distribution of peaks and troughs. In contrast, with the inclusion of ground elevation, 
canopy patterns exhibited anisotropy and had a fractal dimension of 1.78. At the sensor scale, the steepness 


of ground slopes was unrelated to canopy height measures. 
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INTRODUCTION 


Forest canopy structure is comparable to a 
fractal sponge—a volume of spaces and phyto- 
mass between two (as found with a planar ob- 
ject) and three (as found with a solid object) 
dimensions. By ignoring fine resolution spaces, 
our analysis was able to depict the outer surface 
or forest canopy hull. This virtually unbuffered 
surface experiences dramatic changes in micro- 
climate (wind, temperature, humidity) and is in 
part sculpted by erosional forces. The forest can- 
opy hull is dynamic and subject to phenological 
changes (leaf abscission and flowering) and to 
longer term ecological changes (resulting from 
competition and mortality). Its shapes and tex- 
tures embody the history of the forest mosaic. 
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Forming this mosaic are different tree species, 
age classes, underlying soils, ground topogra- 
phy, and disturbances (both endogenous, such as 
treefalls, and exogenous, such as wind throw). 
Though perhaps belittling the complexity of a 
dairy product, forests are not just Swiss cheese, 
as stated by Lieberman et al. (1989); they are 
not merely a random distribution of gaps (or 
holes) of different sizes. The forest canopy sur- 
face can be highly organized into wavelike (Sato 
& Iwasa 1993) or multifractal formations, frac- 
tals embedded within fractals (Solé & Manrubia 
1995, Drake & Weishampel 2000), where the 
gap size distribution follows a power law. Such 
emergent properties, useful in categorizing for- 
est types, also affect canopy-atmosphere inter- 
actions (e.g., energy and mass exchange and 
transport) and ecological processes (e.g., recruit- 
ment, growth, and competition). The unevenness 
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FIGURE 1. 
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25m 


Portion of swath (~1 km?) of 2720 (80 X 34) overlapping LVIS footprints across a region of 


primary forest in La Selva, Costa Rica. The solid line in the inset is a tributary of the Puerto Viejo River. Circles 


represent the distribution of footprints. 


of the top of tropical rain forest canopies may 
promote turbulent eddying of air, thereby en- 
hancing transpiration, photosynthesis, and res- 
piration (Leigh 1999). 

Accurate mapping of forest canopy topogra- 
phy that depicts the spatial distribution of gap, 
dominant, codominant, and emergent crown el- 
evations was deemed impossible 25 years ago 
(Leonard & Federer 1973). Though not impos- 
sible, such mapping remains difficult to obtain. 
Detailed canopy topography for the most part 
has been limited to areas a few ha in size, with 
canopy access provided by construction cranes 
(Parker et al. 1992, Parker 1993, 1995). Laser 
altimeters, originally designed to measure 
ground topography, produced results that were 
confounded by interfering returns from forest 
canopies (Arp et al. 1982). Laser altimeters 
adapted to focus on this signal are now capable 
of mapping canopy structure of forested land- 
scapes (see Weishampel et al. 1996). We ana- 
lyzed patterns of the canopy surface in a ca. 
100-ha section of primary rain forest that were 
derived using the latest forest inventorying laser 
technology. Knowledge of such patterns is in- 
strumental in developing radiative transfer and 
gaseous flux models. 


METHODS 


The Laser Vegetation Imaging Sensor (LVIS), 
see Blair et al. (1999), acquired canopy and 


ground measurements from La Selva Biological 
Research Station in Costa Rica (FIGURE 1), prior 
to the wet season in March 1998. LVIS flew 
aboard the NASA C-130 at altitudes ca. 8 km 
above ground level. LVIS is an imaging laser 
altimeter, similar to its immediate predecessor 
from the NASA Goddard Space Flight Center, 
called the Scanning Lidar Imager of Canopies 
by Echo Recovery (SLICER) system (Blair et 
al. 1994, Lefsky et al. 1999a, 1999b, Means et 
al. 1999). It captures canopy structural patterns 
based on the reflectivity of component surfaces 
(such surfaces as leaves, branches, epiphytes, li- 
anas, vines) from the canopy top to the ground. 
The resulting profile for a canopy is a waveform 
representing the vertical distribution of (approx- 
imately) nadir-projected surfaces. LVIS, unlike 
its predecessors which consisted of narrow tran- 
sects of smailer laser footprints, is capable of 
mapping canopy landscapes (Weishampel et al. 
2000). As programmed for this mission, it pro- 
jected a 1-km across-track swath consisting of 
80 overlapping, 25-m diameter laser footprints 
with a vertical resolution of 30 cm (FIGURE 1). 
Canopy surface height measures differ, de- 
pending on whether height is measured from the 
ground up using traditional forest clinometer 
techniques (Avery & Burkhart 1994) or depth is 
measured to the upper surface of the canopy 
from above using a point drop technique (Miller 
& Lin 1985, Parker et al. 1992). The clinometer 
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FIGURE 2. Contour and surface maps of canopy heights not including ground topography. Grey-scale intervals 


of 8 m range from 4 m (white) to 44 m (black). 


method to measure tree height, typically does 
not account for differences resulting from fluc- 
tuations in ground elevation; it is equivalent to 
measuring the difference between the first and 
last laser return above the background noise 
from LVIS (FIGURE 2). The point drop method, 
analogous to bathymetric techniques, convolves 
ground and canopy elevations; it is equivalent to 
measuring the distance to the first LVIS return 
above the background noise (FIGURE 3). The cli- 
nometer method may be used for estimates of 
aboveground biomass or merchantable timber 
(Nelson 1997); whereas the point drop method 
may be more appropriate for canopy roughness 
measures affecting turbulent air flow (Groß 
1993). Moreover, tree height and ground topog- 
raphy are often related; for instance, trees on 
sloped terrain tend to be shorter than trees on 


level terrain, and trees tend to decrease in size 
as elevation increases. 

To compare topographic patterns of canopy 
height data derived using these two methods on 
the roughly 1 km? area of primary rain forest 
shown in FIGURE 1, we used several techniques. 
First, we followed a hypsometric approach root- 
ed in geomorphology (Strahler 1952). Hypso- 
metric methods, when applied to forest canopies 
(e.g., Leonard & Federer 1973, Lefsky et al. 
1999a, G.G. Parker unpubl. data), depict the 
manner in which the canopy volume beneath the 
upper surface is distributed. These methods pro- 
vide a graphical technique for comparing surface 
structure of canopied landscapes of different siz- 
es and topographic relief. To directly assess the 
differences in the spatial distribution of heights, 
we used a multiple resolution goodness of fit 
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FIGURE 3. 
m range from 60 m (white) to 132 m (black). 


procedure (Costanza 1989). For this procedure, 
heights were divided into quartiles and then 
treated as though they were unrelated categorical 
data. We then evaluated the number of height 
categories in windows of different sizes in cor- 
responding locations from two maps. Autocor- 
relative patterns were calculated with geostatist- 
ical/fractal measures (Palmer 1988) and Moran’s 
I (Legendre & Fortin 1989). Such measures have 
been used to quantify stand-level spatial patterns 
of measured (Ford 1976, Drake & Weishampel 
2000) and simulated (Weishampel & Urban 
1996) tree heights. To test for anisotropy, we 
examined autocorrelative properties from direc- 
tions offset by 90°. Lastly, we evaluated the de- 
gree to which ground slope and tree heights 
were related. Vectors representing differences in 
laser-derived ground elevation were calculated 
(FIGURE 4). These show the direction and mag- 
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Contour and surface maps of canopy heights plus ground topography. Grey-scale intervals of 8 


nitude of slope which could be used to predict 
stream networks. Vector magnitudes at a given 
location were compared to canopy height mea- 
sures. 


RESULTS AND DISCUSSION 


Both measures of canopy heights neared nor- 
mal distribution (FIGURE 5). The mean above- 
ground canopy height for this section of old- 
growth forest was 29.3 m with a standard de- 
viation of 7.2 m (FIGURE 5A). This is ca. 5 m 
taller than measured by Clark et al. (1996) at La 
Selva. Heights ranged 5.1—51.5 m. With the ad- 
dition of elevation above mean sea level (msl), 
the mean height of the outer canopy surface was 
99.1 m with a standard deviation of 9.9 (FIGURE 
5B). These measures ranged 56.5—136.8 m, sim- 
ilar to what would be expected with published 
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FIGURE 4. Vector map of ground topography depicting the direction and steepness of slopes. 


ground elevation (Sanford et al. 1994). Thus ca. 
34 m of canopy surface relief across this area 
could be accounted for by changes in ground 
elevation. 

Both hypsometric curves (FIGURE 6) are sim- 
ilar. They are sigmoidal, indicative of normally 
distributed patterns. Their inflection points, 
which signify the transition point where the rate 
of decrease of canopy volume upwards changes 
from an increasingly rapid rate of decrease to a 
diminishing rate of decrease, occur approxi- 
mately at the same location. The slope at the 
inflection point for the canopy heights without 
elevation is slightly less than that with elevation. 
This reflects the less jagged topography in can- 
opies with elevation, as shown in FIGURES 2 and 
3. This difference is also apparent in compari- 
sons of the fractal dimensions (Zeide 1991, Pa- 
chepsky et al. 1997) of the different measures of 
canopy surface. Canopy topography without el- 


evation had a fractal dimension of 1.96; that 
with elevation was 1.78. These dimensions mean 
that canopy topography without elevation was 
closer to random and occupied more of a two- 
dimensional area. 
Though the multiple resolution goodness of fit 
between the two measures improved as window 
size increased (FIGURE 7), it always was less 
than expected with the random permutation 
maps of the quartile categories. The average fit 
(such as F, from Costanza 1989), based on equal 
weighting of the window sizes, was 0.75 on a 
scale of zero to one, with one representing iden- 
tical maps. This fit was significantly lower than 
that found with random quartile maps (F, = 
0.86). Thus at fine to coarse scales, essentially 
no correspondence was found between the 
height measures across the landscape. 
Autocorrelative patterns for the two topo- 
graphic measures were strikingly different. For 
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FIGURE 5. 


Distribution of canopy height measures derived from the laser altimeter: a) without ground to- 


pography and b) with ground topography. The solid lines are normal distributions from mean and standard 


deviation of measures. 


canopy height measures alone, autocorrelation 
decreased to zero after 100 m (FIGURE 8), sig- 
nifying randomness at these scales. This rapid 
decline and the absence of anisotropy (as des- 
ignated by the similarity between 0° and 90° 
correlograms) are what would be expected from 
canopy height measures at La Selva (Clark et 
al. 1996) and from spatial gap model simula- 
tions near these latitudes (Weishampel & Urban 
1996). For canopy height measures that includ- 
ed ground elevation, autocorrelation was ini- 
tially higher and declined considerably up to 
200 m. In addition, anisotropy was pronounced. 
Ground plus canopy topographic patterns ex- 


hibited periodic undulations at lags >600 m. 
The slight negative relationship found between 
slope steepness and canopy height was not sig- 
nificant. 


CONCLUSIONS 


Because ground topography accounted for 
more than a three-fold increase in average can- 
opy surface height and a nearly two-fold in- 
crease in range, the two measures of canopy 
topography were not expected to be similar. 
Furthermore for this section of forest, ground 
elevation changes are relatively low and 
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FIGURE 6. Hypsometric curves of canopy height 
measures. 


Should only minimally affect forest structure. 
Based on the hypsometric graphs, differences 
appeared to be slight, but they were most ap- 
parent when we examined autocorrelative 
properties. In these analyses, the underlying 
geomorphology revealed autocorrelation at 
distances well beyond that expected with sim- 
ple measures of tree height. Additionally, the 
fractal measures showed canopy surface 
height without ground topography to have a 
rougher texture than when ground topography 
was included. Thus the vantagepoint of a bird 
flying above the canopy differs substantially 
from that of a monkey climbing from tree to 
tree within it. 


Goodness of fit (Fy) 





At certain scales, closed canopy forests re- 
semble a head of broccoli (H.H. Shugart pers. 
comm.); that is, they display self-similar fractal 
properties. Measures of canopy surface topog- 
raphy presently limited to a few areas soon will 
be available globally via remote sensing. The 
Vegetation Canopy Lidar (Dubayah et al. 1997) 
is a Satellite mission scheduled to launch in 
2002. It will provide transects of canopy sur- 
face height and below canopy ground topog- 
raphy data. This characterization of canopy 
properties will be used primarily for terrestrial 
ecosystem and climate modeling. It will pro- 
vide a new means of classifying vegetated sur- 
faces and more accurate parameterization of 
models that require information on canopy 
height. For carbon models that need better es- 
timates of aboveground biomass, canopy height 
measures without topographic ground data will 
be available. For boundary layer airflow mod- 
els that need better estimates of aerodynamic 
roughness or surface heterogeneities (Pielke et 
al. 1998), both measures of canopy surfaces 
may be utilized. Hence, measures of canopy 
structure once thought to be logistically impos- 
sible will be commonplace at landscape scales 
globally. 
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Goodness of fit versus window size for the canopy height surface measures with and without 


ground elevation. The dotted lines are the minimum and maximum values from 100 randomly permuted quartile 


maps. 
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